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ABSTRACT 

We present a model for nonlinear decay of the weak wave in three-dimensional incompressible 
magnetohydrodynamic (MHD) turbulence. We show that the decay rate is different for parallel and 
perpendicular waves. We provide a general formula for arbitrarily directed waves and discuss particular 
limiting cases known in the literature. We test our predictions with direct numerical simulations of 
wave decay in three-dimensional MHD turbulence, and discuss the influence of turbulent damping on 
the development of linear instabilities in the interstellar medium and on other important astrophysical 
processes. 

Subject headings: turbulence, MHD 



1. INTRODUCTION 

The interstellar medium (ISM) is turbulent on scales 
ranging from AUs to kpc (see Armstrong et al. 1995, 
Elmegreen & Scalo 2004) with an embedded magnetic 
field that influences almost all of its properties. Turbu- 
lence is essential for understanding the ISM, the intra- 
cluster medium (ICM), the Earth's magnetosphere, solar 
wind, accretion disks, etc. 

The literature on astrophysical turbulence and its ap- 
plications is vast (see Biskamp 2003, McKee & Ostriker 
2007 and references therein) . The concepts of waves and 
eddies are both used to describe magnetized turbulence. 
It is accepted that weak Alfvenic turbulence consists 
of waves, while for strong turbulence wave-eddy dual- 
ism is important (Goldreich & Sridhar 1995, henceforth 
GS95). In strong Alfvenic turbulence, nonlinear interac- 
tions damp the energy of wave-eddies over approximately 
one wave period. However, there are quite a few sources 
which exist in the turbulent environment that emit waves 
with different degrees of monochromaticity. These waves 
propagate in the turbulent medium while being modified 
by the medium. The evolution of such waves is of con- 
siderable astrophysical interest and has been addressed 
in the literature on numerous occasions. For instance, 
collisional and collisionless damping of fast modes in a 
turbulent medium was shown to differ considerably from 
the case when waves propagate along laminar magnetic 
field lines, with important consequences for cosmic ray 
propagation (Yan & Lazarian 2004). The latter, how- 
ever, is the linear damping of waves, while the focus of 
this paper is the nonlinear interaction of Alfven waves 
with turbulence. 

Since we are addressing the subtle issue of nonlinear 
damping of Alfven waves by surrounding strongly non- 
linearly damped Alfvenic eddies, we feel that a discus- 
sion of the properties of MHD turbulence is due. It 
is well known that magnetohydrodynamics (MHD) de- 
scribes the dynamics of a highly conducting fluid with 
magnetic field. Highly disordered quasi-stochastic flow, 
usually called turbulence, is generally observed in set- 
tings when the Reynolds number, the ratio of the inertial 
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force to the viscous force, is large. When conductivity of 
the fluid is very low, it can be described rather well by the 
Navier-Stokes equations, which are, basically, the MHD 
equations without magnetic field. If the conductivity, 
however, is high, the magnetic field can no longer be ig- 
nored. In this case, the magnetic Reynolds number Rem, 
or the ratio of the fluid velocity to the magnetic diffu- 
sion velocity, is large. In fact, in the turbulent flow of 
a conductive fluid there is a mechanism called dynamo, 
which increases total magnetic energy, making it more 
dynamically important over time (see, e.g., Vishniac et 
al, 2003). It is no wonder that in the astrophysical envi- 
ronment, where both Reynolds and magnetic Reynolds 
numbers are high, magnetic fields are always present and 
always dynamically important. 

The paradoxical nature of turbulence is that even 
though the dissipation parameters (e.g. viscosity and 
magnetic diffusivity) tend to zero, the dissipation rates 
stay pretty much the same. The general explanation for 
this, as proposed by Kolmogorov (1941), is that the en- 
ergy is carried from large scales through many interme- 
diate scales down to small dissipation scales. This prop- 
erty is shared by turbulence in magnetized fluids. Much 
of the progress in understanding turbulence comes from 
so-called Kolmogorov models which assume a cascade of 
a conserved quantity, such as energy, down-scale. Once 
the functional dependence of the cascade decay rate, lo- 
cal in k-space, is known, we roughly know the spectrum 
of the conserved quantity (normally energy, but some- 
times other quantities, such as the number of quasipar- 
ticles, enstrophy, etc. )^. These models do not attempt to 
describe all properties of turbulence, but some important 
ones, such as the spectrum of energy. They essentially 
assume self-similarity, so they don't describe deviations 
from self-similarity, known as intermittency. The level 
of sophistication of these so-called mean field models is 
enough for many applications. 

While the spectrum of hydrodynamic incompressible 
turbulence is well-established, strong MHD turbulence is 

^ The Kolmogorov- type models do not necessarily give —5/3 
spectrum, which is the result specific for the incompressible hydro 
turbulence. For example, so-called Langmuir turbulence has the 
cascade direction inverted and the fiat spectrum of energy. 
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an area of active research (e.g. Montgomery & Turner 
(1981), Shcbalin, Matthacus & Montgomery (1983), Hig- 
don (1984)). Numerical research (Cho & Vishniac 2000, 
Maron & Goldreich 2001, Cho, Lazarian & Vishniac 
2002, Cho & Lazarian, 2002, 2003, Miiller, Biskamp 
& Grappin 2003) is roughly^ consistent with a particu- 
lar model of strong turbulence, called Goldreich-Sridhar 
model (Goldreich & Sridhar, 1995). This mean-field 
model predicts structural anisotropy with respect to the 
magnetic field according to so-called critical balance, 
k±Sv ^ kpjA, as well as a k~^^^ energy spectrum. 

The wave dissipation rate has two major applications. 
First, if it is directly measured, either in experiment, ob- 
servation, or numerical simulation, it can be used as a 
test of turbulence itself, and, in certain mean-field mod- 
els of turbulence, since the latter rely on the phenomeno- 
logically estimated value of the mean dissipation rate for 
every particular scale. 

The other application is turbulent damping of linear in- 
stabilities. If the instability growth rate is smaller than 
the turbulent damping rate, the instability does not de- 
velop. This can have numerous consequences. For in- 
stance, it has been customary to compare the instability 
rate to the viscous or magnetic dissipation rates. Because 
of this, almost any "large scale" perturbation with posi- 
tive imaginary part of the frequency obtained from linear 
analysis has been thought to develop a growing unstable 
mode. It is possible, however, that this instability is slow 
enough to be damped by the turbulent field. This way 
the unstable configuration is stabilized, provided it is not 
disrupted by turbulence itself. In the case when turbu- 
lence itself is generated by an instability, the turbulent 
dissipation rate will provide a powerful generic mecha- 
nism for nonlinear saturation of the instability. Such a 
generic saturation mechanism is highly desirable, since 
typically nonlinear saturation models are very specific to 
the instability in question and contain quite a number of 
assumptions. 

One example of the application of turbulent damping 
rates to astrophysics is the study of the escape of high en- 
ergy Cosmic Rays (CRs) from our Galaxy (Yan & Lazar- 
ian 2002, 2004, Farmer & Goldreich, 2004). Lazarian and 
Beresnyak (2006) developed a model in which the low- 
energy CR mean free path is effectively reduced in the 
presence of compressive motions due to the CR-Alfven 
anisotropic instability. 

Previous works that considered the turbulent wave dis- 
sipation include HoUweg (1984), Similon & Sudan (1989), 
Kleva & Drake (1992), Farmer & Goldreich (2004), Bian 
& Tsiklauri (2007) and others. 

In our paper we develop a model of turbulent dissi- 
pation which is purely nonlinear (does not depend on 
either Re or Rem) and has self-similar power-law depen- 
dencies on the wavenumber, which is characteristic of 
turbulence. We also explain why parallel waves (with 
wavevector parallel to the field) decay in a different way 
than perpendicular waves. We give a general formula 
for the decay rate for arbitrarily directed"^ waves. The 
study of the turbulent decay of parallel directed waves 

^ The particular value of the spectral index and the index of 
scale-dependent anisotropy is still an issue of debate (Galtier, Pou- 
quet & Mangcncy 2005, Boldyrev 2005, 2006, Beresnyak & Lazar- 
ian 2006, Gogoberidze 2007). 

^ We call Alfvenic or pseudo-Alfvenic waves paral- 



by Farmer & Goldreich (2004) was motivated by the fact 
that such waves have the highest instability grows rate 
for CR-plasma instabilities such as the streaming and the 
anisotropy instabilities (e.g., Kulsrud 2005). 

In what follows we describe the model of turbulent dis- 
sipation and derive the formula for the dissipation rate 
in §2, provide numerical evidence for turbulent wave dis- 
sipation in §3, discuss astrophysical implications in §4, 
mention competing mechanisms for wave damping in §5, 
discuss previous work in §6 and summarize our results in 
§7. 

2. WAVE DISSIPATION 
The equations of incompressible ideal MHD, 

dtw+ + {w- ■ V)w+ = -VP, (1) 
dtw- + (w+ • V)w- = -VP, (2) 

V • w+ = V ■ w" = 0, (3) 

written in terms of Elsasser variables w"*" = v + b and 
w~ = V — b, where v is the velocity and b is the mag- 
netic field in velocity units b = B/(47rp)^/^, bear close 
resemblance to the Navier-Stokes equations. This re- 
semblance is somewhat superficial, since w+ and w~ do 
not transform the same way velocity transforms under 
Galilean transformations. Due to this fact, there is usu- 
ally some residual, or local mean magnetic field, that 
cannot be excluded by the choice of frame of reference 
(Kraichnan, 1965). If we go to sufficiently small scales 
the decaying perturbation fields will be much smaller 
than the local mean field. This case is called sub- 
Alfvcnic turbulence. The important conservation laws 
in non-dissipative MHD is the conservation of the inte- 
grals / |w+|^rf^x and / \w~\^(fix, which are the analogs 
of conservation of energy in the Navier-Stokes equations. 
This allows us to build phenomenology of energy trans- 
fer in the spirit of Kolmogorov (1941) for each of the w+ 
and w~ variables with two different dissipation rates'*. 
The setup when those dissipation rates are equal is called 
balanced turbulence, while the opposite case is called im- 
balanced or cross-helical turbulence. 

The major difference between the incompressible 
Navier-Stokes and MHD equations is that the character- 
istics of the latter, due to the presence of local magnetic 
field, are always wave-like, with the speed of the pertur- 
bation being va- While it is not possible to have waves 
in incompressible hydrodynamics, it is possible to have 
waves or wave- like perturbations in incompressible MHD. 
According to the GS95 model, most perturbations pro- 
duced by turbulence cannot be called well-defined waves. 
But there are other perturbations generated, for exam- 
ple, by instabilities, which, as we will demonstrate in this 
paper, have frequency much larger than decay rate and, 
therefore, can be called waves. 

Icl/pcrpendiculajr/arbitrarily directed instead of propagating, 
since these waves always propagate along the local magnetic field, 
regardless of the orientation of their wavevector. 

^ In this paper we consider only the symmetric case in which 
the kinetic viscosity equals magnetic diffusivity (Pr = 1), other- 
wise some exchange of energy between Elsasser fields is possible at 
the dissipation or viscous scale. High Prandtl number turbulence, 
which is often relevant for astrophysical processes, was studied in 
Schekochihin et al 2002, Laaarian, Vishniac & Cho 2004, etc. 
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This paper is focused on the study of the decay of a 
small perturbation in stationary forced turbulence. Un- 
like decaying turbulence, where the energy decays self- 
similarly, or as a power-law with time, the decay of a 
small perturbation in stationary turbulence will be expo- 
nential. This statement is based on the general concept 
of the flow of energy through scales, which underlies all 
Kolmogorov-type models. Indeed, in Kolmogorov mod- 
els, the decay rate is some function of the perturbation 
strength and the length scale: 7turb(<^^', A). For exam- 
ple, in strong hydrodynamical turbulence, 7turb = 5v/\. 
Since the wave amplitude is small, it can not significantly 
affect 5v. Therefore the dissipation rate will be constant, 
meaning that the wave will decay exponentially. 

In this paper we only consider the Alfven and pseudo- 
Alfven waves that exist in incompressible turbulence. 
The study of the decay of all three wave modes in generic 
compressible turbulence is a task that will be addressed 
elsewhere. 

Now we will provide a phenomenological description of 
wave decay and the formula for the decay rate for arbi- 
trarily directed waves. To highlight the general problem 
we first consider parallel waves. According to the second- 
order perturbation theory, such waves are not cascaded 
(e.g., Galtier et al 2002). Does this property necessarily 
apply to strong turbulence? Farmer & Goldreich (2004) 
argued that this is not the case, addressing the issue of 
the wave propagation parallel to the local magnetic field 
when the effective perpendicular wavenumber is obtained 
by so-called field wandering. Let us imagine a plane 
wave with very large perpendicular correlation length, 
corresponding to /c^ = and finite parallel correlation 
length A, corresponding to fcy = 1/A. This wave, mov- 
ing a distance of A, will get disrupted by field wandering, 
as various pieces of this wave propagate along its local 
field direction, which is slightly different along the trans- 
verse spread of the wave. Note that only field- wandering 
from the counter-propagating waves will contribute to 
this process, as the co-propagating waves will move ex- 
actly the same distance, being stationary in the frame 
of the wave. The "effective" perpendicular wavenumber 
fc±.off is determined by the condition fcx.cff = ((5B/B)fc||, 
where we use 5B of counter- wave on scale l/fc_L^efF- 

The arguments above can be generalized for the case of 
an arbitrarily directed wave in the following way. While 
the wave is being decorrelated by turbulence, it obtains a 
new "effective" fcx.cff from the vector sum of the original 
fc_L and the one produced by turbulent field wandering. 
The effective k±^cS in the RMS sense will be the square 
mean of the two. If the original wave has a pitch angle 
6 with respect to the field and a wavenumber magnitude 
of k, its effective perpendicular wavenumber will be de- 
termined by the equation 

fc±,eff = k{sin^ e + {SB/Bf cos^ 6»)^/^ (4) 
or, if we assume the GS95 scaling SB/B ~ kj^^^, 

k±,es = fc(sin2 e + {k^,esL)-y^ cos^ e)V2. (5) 

According to the arguments above, we should include 
only the part of SB /B which corresponds to the counter- 
wave, since co-propagating perturbations do not influ- 
ence the structure of the wave. This is due to the exact 



solution of the incompressible MHD which involves only 
one wave species of arbitrary amplitude. In eq. (5), for 
the sake of simplicity, we omitted the coefficient 1/2 be- 
fore 5B/B which would indicate that only half of the 
perturbations contribute to the decorrelation effect. 

It is easy to verify from eq. (5), that for perpendicular 
modes A:j_,eff = fcsin6' = k (as if field wandering is unim- 
portant), while for parallel modes, fcx.eflf ~ k^^^ (as in 
Farmer & Goldreich 2004). 

The decay rate is determined by the A;_L,eff we calcu- 
lated from (5): 

^ = ^^fc^/3^L-V3. (6) 

Eq. (6) deserves a more detailed explanation. While cas- 
cading of perturbations in weak turbulence have resonant 

character, that is, they require counter-waves with equal 
and oppositely directed k|| (e.g., Ng & Bhattachargee 
1996, Goldreich & Sridhar 1997, Galtier et al 2002), 
strong turbulence involves nonlinear perturbations where 
strict resonance condition is not required. This can be 
demonstrated if we consider a cascaded perturbation as 
a passively advected vector field. This approach is valid 
as long as the amplitude of this sample or test perturba- 
tion is small enough to avoid any back-reaction to fully 
developed, balanced strong turbulence. Now, according 
to critical balance of GS95, the transverse structure of 
this passively advected field is strongly perturbed by 
counter-waves with the same k± after traveling a dis- 
tance of l/fc|| = kJ^^^L^^'^. On the other hand, accord- 
ing to GS95, this is the correlation length of the counter 
wavepacket, so the "next" perturbation will be applied 
to our test field incoherently. Therefore, these pertur- 
bations of the transverse structure of the test field are 
irreversible and can be called cascading in perpendicular 
direction. Our main point here is that the lateral struc- 
ture of the test field is irrelevant for its cascading (imless 
the lateral structure produces the transverse structure by 
field wandering), as long as cascading happens in strong 
GS95 turbulence^. 

The formulae (5) and (6) give the general dependence 
of the decay rate 7 on arbitrarily directed (oblique) waves 
with wavevector k and pitch- angle 6. This dependence 
is, in general, not a power-law of k. There are two lim- 
iting cases of the perpendicular directed and the parallel 
directed waves which give power-laws, namely 

j^ = VAk'/'L-'/', ^^^=VAk^'^L-^'\ (7) 

Fig. 1 shows the dependence of the decay rate 7 on the 
wavevector k for a variety of pitch angles 6 with respect 
to the magnetic field. 

3. NUMERICAL RESULTS 

We have conducted a series of 71 incompressible three- 
dimensional 320'^ MHD simulations with stochastic tur- 
bulent driving. We used a pseudo-spectral code to solve 
Eqs. 1-3 and introduced linear dissipation as an energy 
sink. We used relatively low-order (6th order) hyperdif- 
fusion to extend the inertial range while avoiding a strong 

® This way the test-field approach is diff'erent from the dynamic 
model of GS95 itself, where the lateral structure of the counter- 
waves is essential for back-reaction on co-waves and creating its 
own critical balance. 
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Fig. 1. — The relative decay rate versus wavenumber for different 
angles as obtained from Eqs. 5 and 6. 



bottleneck effect. The turbulent driving supplied energy 
between k = 2 and 3.5 while most of the dissipation hap- 
pened around k = 70. We regarded the interval k = 5 

to 50 as an incrtial interval of turbulence. Simulations 
had an external magnetic field and the driving was tuned 
to produce strong turbulence which was nearly isotropic 
around the driving scale. More details of the code and 
the driving can be found in Cho and Vishniac (2000). 

In this series one of the simulations was without wave 
driving, and all others were driven with a plane wave 
with a particular wave vector k. The turbulent driving 
and the initial conditions in all series were exactly the 
same. The initial conditions were taken from another 
turbulent simulation with driving that ran for around 
20 Alfven times. We waited for about 0.3 Alfven times 
for the wave to reach saturation. We chose the wave 
driving amplitude so that the saturated wave does not 
incur any significant back-reaction to turbulence. After 
this the wave driving was switched off but the turbulent 
driving continued just like in the simulation without wave 
driving. 

As expected, after the wave reached saturation, it was 
no longer plane-parallel and its energy distribution in 
k-space was affected by both field line wandering and 
strong turbulent cascading. It is more appropriate to 
say that the simulation box was filled with a collection 
of wave packets which had the dispersion in the direction 
of k and the dispersion in and k± that comes from tur- 
bulent decorrelation and the corresponding uncertainty 
relation. In Fig. 2 we presented a visual representation of 
the distribution of wave energy that was initially driven 
as the parallel wave with A;|| = 17 (in cube size units) and 
reached saturation. 

The wave energy was obtained by subtracting the re- 
sults of the simulation with and without wave driving. 
We summed the result of subtraction only over the rele- 
vant regions in fc-space, to reduce errors. We note again, 
that driving, although being turbulent-like, was prede- 
termined, and within our short simulation times (about 
0.6 Alfven times) chaos in motions on large scales was not 
yet able to develop. This was confirmed by the fact that 
the result of subtraction was mostly positive and signif- 
icantly differed from zero only around the region where 
we pumped the wave. One notable exception from this 
rule was the velocity at k=0, or the mean velocity. It 
increased, which is interpreted as a transfer of momen- 
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Fig. 2. — Decorrelation of the parallel-propagating wave by tur- 
bulent field wandering. Darker areas indicate larger wave energy. 
The wave is decorrelated and cascaded primarily in the direction 
perpendicular to the magnetic field. 
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Fig. 3. — The evolution of the oblique wave energy during the 
driving and decay stages. The wave driving was at (fc||,fex) = 
(12,12) in code units, approximately (4, 4)L~^. 



tum from the wave to the fluid. We have carried out a 
total of 70 320^ simulations with different values of the 
wave k vector which have covered the inertial range of 
our turbulent simulation from A: = 10 to fc — 40 with 
different angles with respect to the mean magnetic field. 

After the wave driving was switched off, the energy of 
the wave started to decay exponentially (see Fig. 3). We 
measured the decay rate by fitting a decay curve. Fig. 4 
shows the dependence of the decay rates on the k value 
and on the angle of original inclination of wave to the 
magnetic field (the angle 9 between k and B). As we 
see, for waves parallel to the external field the decay rate 
follows a fc^/^ law (see also Farmer & Goldreich 2004). 
The theoretical curves for the different inclinations that 
were fitted into numerical values of 7, unexpectedly, have 
slightly different values for the outer scale L. We believe 
that this could be either due to the fact that the wave- 
lengths that we considered were relatively close to the 
driving scale, or the fact that we simulated trans- Alfvenic 
turbulence. 

4. APPLICATIONS OF TURBULENT DAMPING 

The damping of Alfven waves with k parallel to the 
magnetic field was invoked by Farmer & Goldreich (2004) 
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Fig. 4. — The dependence of the decay rate on the magnitude of 
the wavevector for different angles between k and Bo in a direct 
3D numerical simulation (points), fitted with curves obtained from 
eqs. 5 and 6 (lines). 



in the context of the damping of the streaming instability 
by ambient turbulence (first mentioned in Yan & Lazar- 
ian 2002). They estimated that with the modest level of 
turbulence in our Galaxy, the streaming instability will 
be suppressed for cosmic rays with energies higher than 
about 100 GeV; therefore such cosmic rays can freely 
propagate along the magnetic field lines. Lazarian & 
Beresnyak (2006) have used turbulent damping as one of 
the mechanisms that limits the CR anisotropy instabil- 
ity which occurs in a magnetized fluid in the presence 
of compressible turbulence. This instability greatly de- 
creases the CR's mean free path. However, due to tur- 
bulent damping, there is a high energy limit of around 
1000 GeV (for our Galaxy) for this mechanism. Parti- 
cles with energy higher than this limit will be unaffected 
by the instability and scattered primarily through other 
mechanisms, such as direct interaction with MHD modes 
(Chandran 2000, Yan & Lazarian 2002, 2004, 2007). 

The influence of the turbulent decay to the develop- 
ment of instabilities can lead to the various interesting 
astrophysical effects. Let us consider one example. In a 
recent paper, Everett et al (2007) (henceforth E07) con- 
sidered a model of launching the Galactic wind by CR 
pressure (see also Breitschwerdt et al, 1991). They ar- 
gued that CR pressure is necessary to launch the wind, 
since thermal pressure is too small. One of the critical 
assumptions of this model is that the streaming instabil- 
ity operates effectively and allows CRs to exchange mo- 
mentum with the plasma. Let us now assume that the 
streaming instability is damped by turbulence. Accord- 
ing to Farmer & Goldreich (2004), if the amplitude of 
turbulence is similar to that of our Galaxy, the damping 
of the streaming instability will affect CRs with energies 
higher than 100 GeV, making them effectively discon- 
nected from the background plasma and able to escape. 
If most of the energy (and the momentum) of CRs are 



in the lowest energy CRs, as E07 originally assumed, 
this would present no problem for their model. It is be- 
lieved, that our Galaxy, on average, has a relatively steep 
spectrum of CRs, between the slopes of —2.6 and —3.1, 
justifying this assumption. However, the winds are most 
likely to be launched from the sites of intensive CR pro- 
duction, such as regions with many supernova shocks. 
Some modern theories of shock acceleration with precur- 
sor favor a rather shallow spectrum of the accelerated 
particles, such as shallower than —2 (Malkov & Drury 
2001) ®. CRs with such a spectrum have most of their 
energy and momentum in higher energy particles that 
can freely escape because of damping of the streaming 
instability and, therefore, CR pressure is too weak to 
launch the wind. The existence of the Galactic wind 
would therefore suggest that it is always being launched 
by a CR population with a steep spectrum, or that these 
winds have very low level of turbulence. Since we cannot 
directly observe the spectrum of the CRs near acceler- 
ation sites such information about the slope being shal- 
lower or steeper than —2 can be used to limit theories of 
CR acceleration in shocks. 

The other application of energy dissipation of Alfven 
waves is the problem of coronal heating. It has been 
long known that due to the very small values of vis- 
cous and magnetic dissipation coefficients, ohmic and 
viscous heating are unlikely to dissipate Alfven waves 
significantly. On the other hand, we know that a rela- 
tively small fraction of energy of waves actually escapes 
the coronal region. Turbulent dissipation can deal with 
this problem, since it provides much higher dissipation 
rates than ohmic and viscous heating, and, unlike the lat- 
ter, docs not depend on the Reynolds numbers i?e and 
Rcrn- Parker (1991) noted, however, that in the corona, 
the stochastic component of the field mostly consists of 
waves propagating in one direction, outwards from the 
Sun. In our terminology this is a strongly imbalanccd 
case. While there is certainly a strong imbalance close to 
the photosphere near the source of the waves, the situa- 
tion in the upper solar corona could be closer to balanced 
turbulence due to various mechanisms that allow waves 
to reflect back, such as the geometry of the guiding mag- 
netic field, the parametric instability, etc. Despite coro- 
nal heating still being a poorly understood process there 
are few things that we would like to note with regards to 
our model. First of all, our model does not consider dis- 
sipation in imbalanccd turbulence, due to the fact that 
models of strong imbalanccd turbulence are being de- 
veloped (see Beresnyak & Lazarian 2007 and references 
therein, see Discussion). Secondly, we only considered 
waves with small amplitudes that do not produce strong 
back-reaction to turbulence. We hope, however, that the 
present work can provide insight into a more complicated 
picture that will include imbalanccd turbulence, reflec- 
tion of waves, etc. 

There is also the problem of launching winds, such as 
a stellar wind, by momentum from waves. It is similar 
to the previous problem, in that we need mechanisms 
of wave reflection, in order to generate turbulence, oth- 
erwise waves will freely escape. It is well known that 
the solar wind is strongly imbalanccd close to the Sun, 

^ The average galactic value of —2.6 is then obtained by a variety 
of mechanisms such as re-acceleration, escape, etc. 
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Fig. 5. — Dominating dissipation mechanisms for various wave 
lengths and wave amplitudes in the ICM or hot ISM. Linear coUi- 
sionless (Landau) damping acts only on the pseudo-Alfven mode, 
while nonlinear Landau damping and turbulent damping operate 
on both modes. 



but becomes less imbalanced at larger distances, almost 
reaching cquipartition between waves near the Earth's 
orbit. Roberts et al (1987) suggested that this was due 
to the dominance of the kinetic energy on the outer scale. 

5. COMPETING MECHANISMS OF DAMPING 

While the ohmic and viscous dissipation rates decrease 
very rapidly with increasing wavelength, as A:^, there are 
a number of mechanisms of Alfven wave dissipation that, 
in principle, can compete with turbulent dissipation. The 
dissipation mechanisms can be subdivided into linear and 
nonlinear mechanisms. We call it a linear mechanism 
when the decay rate docs not depend on the amplitude 
of the wave. Such mechanisms include viscous damping, 
turbulent damping (even though the turbulence itself is a 
nonlinear process), collisionless (Landau) damping, and 
the resonance damping. Nonlinear damping mechanisms 
have a decay rate that depends on the wave amplitude. 
These include nonlinear Landau damping, wave steepen- 
ing, etc. 

Nonlinear Landau damping (e.g., Kulsrud 2005) re- 
quires two waves with close frequencies and a population 
of ions that is in resonance with the beat wave of those 
two waves. While, generically, the decrement of this 
damping is proportional to the frequency, it has a nonlin- 
ear {SB/Bof factor that limits damping to high ampli- 
tude waves. Note that the turbulent damping mechanism 
that we consider in this paper has a substantial limita- 
tion in that the wave amplitude has to be small enough to 
not back-react on the turbulence itself. Thus, nonlinear 
damping mechanisms and the linear (with respect to the 
wave amplitude) turbulent damping are, in a sense, com- 
plimentary to each other (see fig. 5). Another nonlinear 
mechanism of damping is wave steepening (e.g., Cohen 
& Kulsrud, 1974). This mechanism is also supposed to 
have 7 ~ kvAiSB/Bo)"^. 

The prominent linear mechanism of damping is so- 
called resonance damping (e.g., HoUweg, 1984). Unlike 
previously mentioned mechanisms that are considered in 
a homogeneous medium and usually work equally well in 
any geometry, this linear damping mechanism appears 
when, due to the specific boundary conditions or inhome- 
geneity of the flow, the so-called Alfvenic resonance ap- 
pears. This resonance corresponds to the flow speed be- 



ing equal to the phase speed of the Alfvenic wave. The 
resonance is typically the l/x type divergence for the 
magnetic field and velocity. Notably, the dissipation pro- 
vided by this resonance does not depend on viscosity or 
magnetic diffusivity (Kappraff & Tataronis 1977). The 
dissipation rate can be calculated from the ideal MHD 
equations with v = rj = using contour integrals (e.g., 
Livshitz & Pitaevskii 1981, Pariev & Istomin 1996). The 
Alfven resonance of the flow has been proposed as a 
prominent candidate for an in-situ acceleration mecha- 
nism in jets (Beresnyak, Istomin fc Pariev 2003). Reso- 
nance damping usually provides damping rates that are 
a fraction of the wave frequency. It is considered to be 
important when the wave length is of order the scale at 
which boundary conditions are set. In jets we expect 
turbulent damping to be more important than resonance 
damping (if the latter is present) for wavelengths that 
are much smaller than the jet's diameter, assuming that 
for these wavelengths the medium can be treated as ho- 
mogeneous. 

6. DISCUSSION 

This paper deals with wave dissipation in a balanced 
MHD turbulence. The more general imbalanced case was 
not considered for a variety of reasons. First, the the- 
ory of strong imbalanced MHD turbulence is still being 
developed (Lithwick, Goldreich & Sridhar 2007, Beres- 
nyak & Lazarian 2007). On the other hand, if we adopt 
the model of imbalanced turbulence presented in Beres- 
nyak & Lazarian 2007 (furthermore BL07) a number of 
questions arise. Most importantly, we will be unable to 
reproduce the statement of §2 that the cascading of the 
wave is irrelevant to its lateral structure. This is due to 
the fact that neither the strong nor the weak wave from 
BL07 has a critical balance with itself (i.e. k^^VA = k±Sv 
does not hold if all quantities refer to the same wave, in- 
stead BL07 provides a different expression for the critical 
balance of each wave). Therefore there is no guarantee 
that we can use the strong cascading formula, at least 
when the wave has a large A;|| . 

An attempt to include irregularities of magnetic field 
lines due to turbulence to increase predicted dissipation 
in solar corona was made in Silimon & Sudan (1989). 
Their model relies on a description of turbulent mag- 
netic field lines as divergent with a specific Lyapunov 
constant, or exponentiation rate A. We argue that this 
is not a proper description of turbulent fields in devel- 
oped turbulence because such fields are approximately 
self-similar and do not have any designated scale in its in- 
ertial interval. In developed turbulence every scale has its 
own exponentiation rate. Silimon & Sudan (1989) give a 
formula for dissipation length which is inversely propor- 
tional to the Lyapunov constant (which is not specified, 
but estimated) and only weakly (logarithmically) propor- 
tional to the wavenumber. It also depends on magnetic 
Reynolds number. This paper, in contrast, advocates 
purely nonlinear dissipation rates that do not depend 
on Rcm or Re. Parker (1991) criticized Silimon & Su- 
dan's approach of using stochastic field lines to increase^ 
dissipation by noting that the stochastic component of 
the field is itself part of the waves propagating outwards 
along the field. We commented on this controversy pre- 
viously by noting that one has to have a mechanism to 
reflect the waves back, otherwise they escape freely. 
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Kleva & Drake (1992) have studied nonlinear dissipa- 
tion of waves in the presence of a stochastic field on an 
outer scale by numerical methods. They confirmed that 
the dissipation of large scale waves does not depend on 
dissipation coefficients, while for small scale waves it ap- 
proached the viscous and resistive limit. However, their 
stochastic field was not turbulent, but rather a predeter- 
mined field on the outer scale. Also, the damping rate 
they measured did not depend on wavevector as a power- 
law. 

Bian & Tsiklauri (2007) have considered mixing of 
Alfvcn wavepackets in chaotic magnetic fields. They 
obtained an analytical expression for the advection- 
diffusion of wavepackets using a WKB approximation 
and claimed that the stationary wave energy spectrum 
is k~^. We feel that the WKB method is not the appro- 
priate tool to derive turbulent dissipation, as it requires 
that the wave should have wavelengths that are much 
smaller than the underlying perturbations, while in tur- 
bulence the most effective nonlinear dissipation comes 
from scales comparable with the wavelength. 

Hollweg (1984) considered dissipation of Alfven waves 
in coronal loops, primarily due to resonance mechanism, 
but also speculated on turbulent heating. He took the 
Kolmogorov dissipation rate (which is a decay that hap- 
pens on a kinetic timescale X/v) with the perturbation 
length scale as wavelength and the RMS wave velocity 
as perturbation velocity. This way he obtained dissipa- 
tion which is compatible with heating in coronal loops 
as well as coronal holes. We have to note, however, that 
propagating waves do not necessarily generate turbulence 
and even if there is a flux of waves in both directions, 
there are certain requirements for MHD turbulence to 
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be strong turbulence and have the Kolmogorov dissipa- 
tion rate. For example, in the context of coronal loops, 
with a strong guide field, it is easy to imagine a situa- 
tion where turbulence is not strong on the outer scale, 
therefore it has dissipation rate lower than Kolmogorov. 

7. SUMMARY 

In this paper we demonstrated the following: 

1. Parallel and perpendicular directed waves are dis- 
sipated in different ways. Perpendicular waves are cas- 
caded naturally by counter-waves, while parallel waves 
are first decorrelated by the field wandering and then 
cascaded. 

2. In the inertial interval of sub-Alfvenic turbulence, 

parallel waves are damped more slowly than perpendic- 
ular waves with the same wavenumber. The decay rate 
for the parallel wave is much smaller than its frequency, 
it is not eddy but a well-defined wave. 

3. In a homogeneous medium with open boundaries 
turbulent dissipation is more effective than the other dis- 
sipation mechanisms, provided that the wave amplitude 
is sufficiently small and the wavelength is in the inertial 
interval of turbulence. 
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